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Abstract 

We perform Brownian dynamics simulations for studying the self-diffusion in two-dimensional (2D) 
dusty plasma liquids, in terms of both mean-square displacement and velocity autocorrelation function 
(VAF). Super-diffusion of charged dust particles has been observed to be most significant at infinitely small 
damping rate 7 for intermediate coupling strength, where the long-time asymptotic behavior of VAF is 
found to be the product of and exp {—^t). The former represents the prediction of early theories in 
2D simple liquids and the latter the VAF of a free Brownian particle. This leads to a smooth transition 
from super-diffusion to normal diffusion, and then to sub-diffusion with an increase of the damping rate. 
These results well explain the seemingly contradictory scattered in recent classical moleculai^ dynamics 
simulations and experiments of dusty plasmas. 

PACS numbers: 52.27.Lw, 52.27.Gr, 66.10.cg 
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Diffusion processes play an important role in determining dynamical properties of many phys- 
ical, chemical and biological systems UJ]. Diffusion in two-dimensional (2D Vphysical systems is 
of particular interest because of its significance in fundamental physics flQ^Q], as well as its ap- 
plications in many 2D systems, such as strongly-correlated electrons on the liquid helium surface 
||5|], colloidal suspensions core-softened fluid [|7D, strongly coupled non-neutral plasmas [|8D 



and, in particular, strongly coupled dusty plasmas (SCDPs) 190, which have been generally recog- 
nized as a promising model system to study many phenomena in solids, liquids and other strongly 
coupled Coulomb systems at the kinetic level and have been used frequently to study many physi- 



cal processes in 2D systems^ such as heat condu ction [llOh. shear flow and certainly diffusion 



of charged dust particles yJi [13|, ll4ill5Lll6L 
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Diffusion is usually described by using particles' mean-square displacement (MSD) and ve- 
locity autocorrelation function (VAF), defined by MSD(t) = {\ri(t) - ri(O)p) and Z(t) = 
(vj(t) ■ Vj(0))/(t'j^(0)), respectively. Here ri{t) and Vj(t) are the position and velocity of ith 
particle at time t, respectively, and the angular bracket denotes an ensemble average. Note that 
the (normal) diffusion coefficient D is determined by the MSD through the Einstein relation 
D = limi^ooMSD(t)/(4t) and by VAF through Green-Kubo formula D = Z{t)dt/{2T) 
111]. What makes things interesting in 2D systems probably lies in the fact that the diffusion 
processes could exhibit some anomalous behaviors. For the normal diffusion, one would expect 
MSD(t) oc t, and \Z{t)\ decays faster than t~^, when t ^ oo, in order to obtain a meaningful 
diffusion coefficient. Diffusion that do not comply with these characteristics are usually regarded 
as anomalous one, for example the so-called super-diffusion and sub-diffusion. In the former case 
(super-diffusion), one expects 121 bi |4|] MSD(t) oc t^^°', where a is positive, and \Z(t)\ oc 
when t ^ oo and in the latter (sub-diffusion), a becomes negative, when t ^ oo. 

Early molecular dynamics (MD) simulations with 2D hard-disk liquids [21] observed an asymp- 
totic decay of the VAF in the form that is oc which led to the conclusion that particle motions 
in 2D hard-disk systems are non-diffusive (or super-diffusive). Later theories of 2D simple liquids 
I13i|4i] came to the same conclusion, but with a wider argument saying that this asymptotic behavior 
should also exist in kinetic parts of correlation functions for the shear viscosity and heat conduc- 
tivity, regardless of the liquid density and interaction force, and finally lead to the breakdown of 
the hydrodynamics in 2D systems. 

The anomalous diffusion of dust particles in quasi-2D and 2D dusty plasma liquids (DPLs) has 
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been frequently reported in many laboratory experiments 1121 llJ, 114], [15L [16D. In particular, the 
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super-diffusion of char ged dust grains in 2D (monola yer) experiments has been recently observed 



by different groups 1113 



Quinn and Goree [|l3l l were the first to report the observation of 



super-diffusion at long timescales in a 2D DPL experiment, followed by Ratynskaia et al. 1141 
and more recently by Liu and Goree [[^.[n.fisll. These observations were supported by the early 
theoretical pred iction 12, |3|, |4|] and recent MD simulations with 2D and quasi-2D Yukawa liquids 
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However, seemingly contradictory observations exist in both experiments and simulations. 
Nunomura et al. investigated the self-diffusion in 2D DPLs for a wide range of the Coulomb 
coupling strength, and observed a normal diffusion process till the melting point, at which charged 
dust particles are frozen and therefore a sub-diffusion was observed. Their observation was sup- 



ported by their own MD simulation results, and recent Brownian dynamics simulation [|22|], in 
which the neutral damping effect was considered and normal diffusion was observed in the whole 
liquid region. Yet, 2D simulations with softer interactions [|7|] observed normal diffusive motion. 

The present work is particularly motivated by these contrasting observations in recent 2D dusty 
plasma experiments and simulations, and is intended to clarify some of the above incompatibility, 
for example, the effects of the neutral gas damping and stiffness of the dust particle interaction on 
super-diffusion, by using extensive numerical simulations. 

Our simulation is based on the Brownian dynamics (BD) method [23|, which brings the simu- 
lation closer to real dusty plasma experiments, as the effects of the neutral gas damping and dust 
particle Brownian motions are included in a self-consistent manner. Most relevant experiments are 



m. 



those conducted in equilibrium conditions |[l3, 14, 

In our simulation, = 20000 particles are placed in a square with periodical boundary con- 
ditions (PBC). A 5t/i-order Gear-like Predictor-Corrector algorithm [|24ll was used to integrate the 
Langevin equations of every interacting Brownian dust particles. As usual, we assume here pair- 
wise Yukawa interaction between charged dust particles: $(r) = (Q'^/r) exp (— r/A/)), where Q 
is the charge on a dust particle, r is the interparticle distance, and Ad is the Debye screening length. 
Such a system can be fully characterized by three parameters [251]: the Coulomb coupling parame- 
ter F = Q'^/(aT), the screening parameter k, = a/ A, and the damping rate 7 / cu^d due to the neutral 
gas, where T is the system temperature (in energy units), a = (vrn)^^/^ the Wigner-Seitz radius, 
and n the equilibrium dust number density. The dusty plasma frequency is = [2Q'^/ (ma^)] 
where m is the dust particle mass. The system size is approximately 250a x 250a, which together 
with PBC and the dust-acoustic speed [|26|] decides our maximum independent observation time to 
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be 320cj-|. 




FIG. 1: MSD(t)/t (normalized by Wpja^) for k = 1.0, = 0.001 and different T. Dash-lines are 

linear fits (in log-log scale) of asymptotic behaviors at long time limit, a is the slope of the fits. Lengths of 
the dash-lines show roughly the fitting range. Typical uncertainty is ±8% of the shown values. Circles ai^e 
corresponding diffusion coefficients obtained through Green-Kubo formula. 

Our analysis will be based on MSD and VAF. For a better visualization of different types 
of diffusion processes, we have used throughout this paper the quantity MSD(t)/t, rather than 
MSD(t) itself. Let us firstly see the evolution of MSD in the zero damping limit. Figure [T] displays 
MSD(t)/t for K = 1.0 and 'j/uJpd = 0.001, but for different F, which covers from nonideal 



gaseous state (F = 1) to the melting point (F = 180) 112611 . Shown together in this figure are linear 
fits (in log-log scale) of their asymptotic behaviors at short and long time limits, respectively. 
The former gives a constant value of a = 1, indicating the ballistic motion, whereas the latter 
varies for different F, as depicted by the values above or below the corresponding dash lines. 
The variation of a with F at long time limit shows clearly a hill shape, with a flat peak of a 
in the range of F G [50,80] and a peak value of a ^ 0.18. This behavior indicates a clear 
transition from normal diffusion — super-diffusion normal diffusion, when F goes from 1 
to 180. Super-diffusion is most significant at the intermediate Coulomb coupling strength. The 
general tendency agrees qualitatively with previous MD simulations of Ott et al. (in the 2D limit) 




FIG. 2: The velocity autocon-elation function for k = 1.0, 7/wpd = 0.001 and different P. The decay of 
and are also shown for comparison. Dash-line follows t^^ exp {—jt). 



UM and Liu and Goree [|18l] . except that in the latter the peak appears around F = 20 with a value 
of a ~ 0.3. This discrepancy could have been caused by the limited observation timescale (which 
then was determined by the finite size of the system.) in Ref. |18]. Their latest simulations [ 1^ 
with substantially more particles gave a peak value very close to ours. In addition, the diffusion 
coefficients calculated by using the Green-Kubo formula were also shown as circles for T with a 
normal diffusive behavior. Their horizontal location indicates the integration limits. We found a 
very close agreement on the diffusion coefficient by using these two methods, which proves the 
consistence of our simulation. 

Figure [2l depicts the VAFs (their absolutes) for different F in the same conditions as Fig. [U 
Shown together in this figure are different decay rates following t^^-^, t^^ and exp (— 7t), 
respectively. One observes a continuous transition from a fast decay at F = 1 (faster than t~^^) 
to a slow decay at F = 50 approximately following t^^, indicating the occurrence of the super- 
diffusion. The last three tails all follow roughly t^^ decay, however, they seem to be better fitted 
by t^^ exp (— 7t), as shown by dash lines in the figure, with the exponential part coming from 
the VAF of a free Brownian dust particle. Since the exponential factor obviously suppresses the 
super-diffusion, what we observe here is a super-diffusion in a transient timescale. In other words. 




FIG. 3: The velocity autocorrelation function for k = 1.0, F = 10, and different damping rate 7 (normalized 
by Wpd). Heavy dash-lines follow exp (— 7t). The decay of is also shown for comparison. 

there should be no absolute super-diffusion (in a sense of infinitely long time limit) in a damped 
system in the equilibrium state. However, for this specific damping rate = 0.001, its effect 

is negligibly small at our maximum timescale. One would need a timescale of order 10^^^"^^, 
which then require about one million dust particles, to observe a complete suppression of the 
super-diffusion in this case. Nevertheless, this tendency will be clearer next when we vary the 
damping rate. 

Figure [3] displays the VAF for k = 1.0 and F = 10 but different 7, as a typical example 
exhibiting the damping effect on the long time tails of the VAF. The long time decay of all VAFs 
can be well fitted in the form of exp (— 7t) up to about 7/u^pd = 0.1. For higher damping 
rates, the VAF reaches the noise level (about 10^^) too fast, so that no reliable fits can be applied. 
It should also be noted that a damping as small as 7/c<jpd = 0.01 is sufficient to bring the decay 
close to t~^'^ in our timescale and to suppress the super-diffusion. Since the typical value in 
similar dusty plasma experiments [14, l^lll ] is about a few percent of iOp^, we would expect that 
there should be no super-diffusion in the equilibrium state. This explains well the absence of the 
super-diffusion in the experiment of Nunomura et al. [12 11]. 

Figure [4] exhibits the variation of MSD(t)/t for k = 1.0 but for different F and 7 values. When 
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FIG. 4: MSD(t)/t (normalized by Wpja^) for k = 1.0, but different F and damping rates (normalized by 
Wpd). Circles ai^e con^esponding to diffusion coefficients obtained through the Green-Kubo formula. 



7 goes from infinitely small to a finite value, another transition from the super-diffusion — > normal 
diffusion sub-diffusion is clearly seen at long timescale especially from panels with higher F, 
which implies an enhancement of the caging effect due to the neutral damping. This explains the 
sub-diffusion close to the melting point observed in the experiment of Nunomura et al. bill . 

Circles in Fig. |4] correspond to the diffusion coefficients D obtained through the Green-Kubo 
formula, which agree with those from MSD (last points of MSD(t)/t) very well. Relative big 
discrepancies are found for the highest Coulomb coupling strength F = 150 and highest damping 
rates 'j/uJpd = 0.5 and 1.0, which implies that our timescale is not long enough to reach the steady 
value. 

One sees from both circles and long time behavior of MSD(t)/t that D decreases with the 
increase of the damping rate. This is physically clear because the damping effect retards the 
dynamics of charged dust particles. This tendency also agrees qualitatively with observations 



in 3D BD simulations U2, 



27|] and in experiments [|25l] of dusty plasma liquids, in which D was 
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FIG. 5: MSD(f)/t (normalized by WpdO^) for Tcff = f{K)T = 50 (/(k) is given in Ref. IM), j/copd 
0.001 and different k. The inserted figure shows exponential a obtained by fitting at long time limit for 
different T^s, with other conditions the same. The uncertainties in the fits are ±13% of the shown values 
for Tcff = 10 and ±8% for T^s = 50 and 100, respectively. 



found to decay as l/jl + 'j /tOpd) with the damping rate, in a certain range of the Coulomb coupling 



parameter [1221 . 
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Finally, we have investigated the effect of stiffness of dust particle interactions on the diffusion 
processes. Here we vary the screening parameter k, while keeping the effective Coulomb coupling 
parameter Fcfr = /(/«)r constant, where the universal scaling function /(k) is obtained by fitting 
the MD simulation results [|26n . Figure [5] shows the variations of MSD for Fcs = 50, 'j/cOpd = 
0.001, but for different k values, as examples. It is seen that the dust particle motion at this 
coupling strength (and timescale, to be more precisely) is super-diffusive. When k, increases, the 
interaction becomes stiffer and stiffer and more and more like a hard-disk system. During this 
course, one observes that the migration of dust particles becomes weaker and weaker due to the 
decrease of the temperature, except at the highest k, where the diffusion is enhanced shortly after 
the ballistic regime. The effect on the super-diffusion seems complicated, as is shown (the inserted 
figure) by the variation of exponential a at long time limit for different Fcff values. Varying k has 
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little effect for Fcfr = 10, as the a rises slightly from 0.13 at k = 0.5 to 0.138 at k = 4.0. For 
Feflf = 50 and 100, the largest value of a was observed at highest k, indicating the most significant 
super-diffusion. The observation seems contradictory with the theoretical prediction of Ernst et al. 
yD, in which they predicted that the super-diffusion motion should exist in all 2D liquids regardless 
of the interaction force, but in accordance with a previous speculation based on many computer 
simulations that the super-diffusion might be more significant in systems with stiffer interactions. 
However, we did not observe a monotonous increase of a with an increase of k (hence stiffness 
of interactions), which indicates that dependence of the super-diffusion on the interaction force 
might be much more complicated than one would have thought. 

In summary, we have studied the self-diffusion processes in 2D dusty plasma liquids by us- 
ing the Brownian dynamics simulation. The latter shows, for the first time, the transition from 
super-diffusion to normal diffusion then to sub-diffusion with the increase of the damping rate. 
In particular, our observations suggest that there should be no real super-diffusion at long time 
limit in a 2D dusty plasma liquids in the equilibrium state, such as the experimental observations 



by Nunomura et al. piQ. Most of other experiments of super-diffusion were made in systems 
that were not in a strict equilibrium state and therefore collective excitations generally cannot be 
ruled out. This applies to the experiments of Ratynskaia et al. lll4l ll5n. where super-diffusion was 



observed in a viscoelastic vortical dust fluid, and to the experiments of Liu and Goree Hid 
which were made in a driven dissipative system. Our simulation results, therefore, may stimulate 
further experimental investigations. 
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